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The exploration of exotic shapes and properties of atomic nuclei, e.g., œ cluster and toroidal shape, is a 
fascinating field in nuclear physics. To study the decay of these nuclei, a novel detector aimed at detecting 
multiple alpha-particle events was designed and constructed. The detector comprises two layers of double-sided 
silicon strip detectors (DSSD) and a cesium iodide scintillator array coupled with silicon photomultipliers array 
as light sensors, which has the advantages of their small size, fast response, and large dynamic range. DSSDs 
couple with cesium iodide crystal arrays are used to distinguish multiple alpha hits. The detector array has a 
compact and integrated design that can be adapted to different experimental conditions. The detector array was 
simulated using Geant4, and the excitation energy spectra of some alpha-clustering nuclei were reconstructed 
to demonstrate the performance. The simulation results show that the detector array has excellent angular 
and energy resolutions, enabling effective reconstruction of the nuclear excited state by multiple alpha particle 
events. This detector offers a new and powerful tool for nuclear physics experiments and has the potential to 
discover interesting physical phenomena related to exotic nuclear structures and their decay mechanisms. 
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I. INTRODUCTION 


The study of exotic shapes of nuclei is an exciting field that has opened new avenues of research in nuclear physics. Un- 
derstanding the shapes of nuclei is crucial for gaining insight into the underlying physical phenomena governing the nuclear 
structure and behavior. When nuclei are excited to high energies or exhibit a high angular momentum, they can exhibit a variety 
of exotic shapes, such as linear-chain, toroid, cylinders, and bubbles, which are not emerge under normal conditions. The 
toroidal structure of the nucleus was first proposed by Wheeler [1]. Wong systematically studied the conditions for the exis- 
tence of toroidal nuclei in medium-and heavy-mass regions [2—4] and light-mass regions [5, 6]. Recently, various sophisticated 
Hartree-Fock (HF) microscopic methods have addressed the issue of light toroidal nuclei [7, 8]. These theoretical studies sug- 
gest that these exotic shapes arise because of the interplay between nuclear, centrifugal, and Coulomb forces. Exotic resonance 
peaks at very high predicted excitation energies in the 7a disassembly of ?8Si were recently observed, which matched well 
with the excitation energy of the toroidal nuclei [9-12], indicating the successful population and detection of toroidal high-spin 
isomers. 

The Hoyle state has important implications for the nuclear reactions and nucleosynthesis processes that occur in stellar 
environments. Be and !?C nuclei exhibit distinct cluster structures, such as the well-known Hoyle state of 12C and the 
a+?n-+?n cluster structure of He, which closely resemble the 3a-condensate-like structure of the Hoyle state [13]. The Hoyle 
state is the key to understanding the nature of nuclear forces and nuclear structures. Besides, some studies have demonstrated 
the emergence of a 7-bond linear-chain molecular rotational band in 14C [14]. Recently, new evidence for a predicted Hoyle- 
like structure in 160 was found [15]. Research on multiple clustering configurations in 74Mg yielded significant results [16]. 
In addition, typical clustering structures in !*Be and linear-chain clustering structures in neutron-rich 16C have been observed 
[17, 18], indicating that clustering is a general phenomenon observed in light nuclei [19-22]. Therefore, it is plausible to assume 
that heavier conjugate nuclei have similar cluster states. 

The excited states of clustering and toroidal nuclei tend to form exotic shapes that can decay into multiple alpha clusters. 
Experimental investigation of nuclei with alpha-cluster states requires the precise measurement of multiple alpha particles 
emitted during the decay process [23, 24]. To study these phenomena, advanced experimental techniques have been developed 
that can fit the requiement of the measurement of multiple alpha-particle events with high resolution. 


* This work is supported by the Strategic Priority Research Program of Chinese Academy of Sciences (No. XDB34030000), the National Key Research and 
Development Program of China (No. 2022YFA1602404), National Natural Science Foundation (No. U1832129, No. 11975210), and Youth Innovation 
Promotion Association CAS (No. 2017309). 

+ Corresponding author, CAO xiguang, Email:caoxg @sari.ac.cn 


The coincident detection of these particles reveals crucial insights into the fragmentation process of nuclei, which helps in 
understanding the underlying physics of the nuclear structure and decay mechanism. Several detectors have been specifically 
developed for measuring multiple alpha waves and have been used in various research facilities worldwide. One method of 
studying nuclear reactions involving multiple alpha particles is to use detectors that can measure the energy and direction of 
each alpha particle. The Cylindrical Array for Tracking and Spectroscopy (CATS) in the Grand Accélérateur National d’Ions 
Lourds (GANIL) consists of two low-pressure multiwire proportional chambers that can detect and identify alpha particles at a 
high counting rate [25]. A state-of-the-art 47 array of charged-particle detectors called ChAKRA works at the Variable Energy 
Cyclotron Center, which facilitates high-resolution charged-particle reactions and spectroscopy studies [26]. The FAZIA can 
detect the Fermi energy domain of charged particles based on three telescope stages [27]. CSHINE, a detector for studying 
the state of asymmetrical nuclear matter, can offer opportunities for experimental studies on the collision dynamics and nuclear 
equations of state in heavy-ion reactions at Fermi energies [28]. Peking University teams performed calibration tests on two 
annular detector arrays and produced outstanding results in discriminating light-charged particles [29]. 

The measurement of multiple a particles is a hot topic in nuclear physics because it can provide new insights into the exotic 
structure and dynamics of the nuclei. Multiple a particles from clustering nuclei decay are challenging to measure owing 
to the limitations of the solid-angle covering and low reaction cross-section. High-resolution and high-sensitivity detectors 
are required to reconstruct the breakup process of nuclei that emit multiple a-particles. However, these detectors are rare and 
expensive, and they must deal with a large amount of background information that can obscure rare events of interest. Therefore, 
new methods are required to reduce background interference and enhance the signal-to-noise ratio [30-32]. The development of 
specialized detectors and innovative data analysis techniques is essential for advancing this field of research. Recently, machine 
learning methods have been used to study and analyze Hoyle states [33, 34] and clustering structures [35], which can be used 
to classify and predict experimental data. The future of this field lies in the design of compact and dedicated detectors that can 
address specific experimental difficulties and challenges, particularly in sophisticated experiments used to explore interesting 
clustering phenomena in nuclei [36-39]. The following sections introduce the specific structure and detection performance of 
the detector in terms of design, construction, simulation, and measurement. 


Il. DESIGN AND CONSTRUCTION 


The detector consists of two layers of double-sided silicon strip detectors (DSSDs) and CsI detectors with a frame arranged 
to maximize the sensitive area for particle collection and detection. Frame design is crucial, as shown in Fig. 1. Telescope 
arrays are built with aluminum frames because aluminum is low-cost, lightweight, and the reaction and activation cross sections 
of aluminum are low; therefore, the impact on detection is small. Therefore, the front frame of the cesium iodide is made of 
high-purity aluminum. 

The overall detector structure can be broadly divided into four parts. As shown in Fig. 1, the gap between the two DSSDs of 
type BB7 is 10 mm. To prevent the DSSDs from being damaged by the tension of the Kapton cables, four aluminum pads are 
attached to secure the Kapton cables. 

A custom-made frame fits the trapezoidal CsI tightly and protects it from damage. The frame also contains a slot for the 
signal wire of the DSSDs. The cesium iodide crystal is sandwiched between a reflective layer and an aluminized film on the 
front surface to collect light. 

The optoelectronic conversion readout model consists of 25 silicon photomultipliers (SiPMs) and their power supply and 
control circuits. Each SiPM is attached to a CsI crystal. The frames ensure that the SiPM and CsI are well aligned so that the 
scintillation from the CsI can be collected and converted into a charge signal by the SiPM as much as possible. The frames 
also provide mechanical support and stability to the module and prevent the CsI from sliding out of position. The SiPM signal 
readout model incorporates a temperature feedback circuit that can correct the signal amplitude due to the temperature shift. 

The signals from the four 34-channel signal cables from DSSDs and 25-channel signal lines from SiPMs are collected using 
a circuit board connected to the electronics analysis module. The 25-channel signal lines are split into common 16-channel 
signal cables that can match various data acquisition modules. 


A. Detectors 


In particular, a Si-CsI combination detector for multiple a particle coincidence measurements can detect the energy and 
position of multiple œ particles with high precision [40—42]. This can help to explore the fundamental laws of nature governing 
the interactions of these particles. This design utilizes 25 CsI crystals, each of which is a sloping ridge with a hypotenuse 
pointing towards the collision center, causing the a particles to deposit energy in only one crystal. The front and rear sides of 
each crystal measured 12.96 x 12.96 mm? and 15.16 x 15.16 mm?, respectively, at a height of 50 mm. The 5 x 5 CsI crystals 
were sandwiched between TeOz for good reflection and combined into the frustum of a square pyramid with dimensions of 
66 x 66 mm? on the front side and 77 x 77 mm? on the rear side. An array of 25 CsI crystals was machined from a large 


Fig. 1. Exploded view of the charged particle detector and the view of the telescope array 


piece of crystal to ensure that their physical properties were similar, and TiO2 was sandwiched between the crystals to act as a 
reflection layer and glue them together. This design renders the telescope array robust and compact. 

The detector consists of three layers: two DSSDs and a CsI scintillator. The BB7-type DSSD from Micron Semiconductor 
Ltd was used to measure the position and energy loss of the incident particles. BB7 type DSSD, which has a larger sensitive area 
and narrower strips than the W1 type, has an active area of 63.96 x 63.96 mm, which consists of 32 strips, approximately 2 mm 
wide, on each side. DSSDs with thicknesses of 300 mm, 500 mm, and 1000 mm were purchased depending on the experimental 
requirements. This type of detector is suitable for high-resolution measurement of charged particles. Heavy ions deposit full 
energy in the DSSDs, whereas lighter particles such as alpha punch through the DSSDs and hit the CsI, where the residual 
energy is deposited. This configuration allows the detector to identify and locate multiple alpha particles simultaneously, which 
is suitable for various applications that require a high angular resolution for charge particle detection. 

The semiconductor material of the Si detector forms a p-n junction, which makes it highly sensitive to charged particles. 
The detector can also locate the position of the particle interaction within the detector with a high spatial resolution, which 
is essential for accurately distinguishing multiple hit events. Moreover, the detector has a low threshold level that allows the 
detection of very low-energy particles. 

By combining the signals from both layers of DSSD and applying kinematic judgments, the particle momentum can be 


accurately determined. The hit position can also be used in conjunction with the final CsI scintillator array to reduce uncertainty 
in the direction of the particle momentum. This can significantly improve the accuracy of particle momentum determination by 
reducing the effects of multiple hits and detector resolution. 

A DSSD with two layers allowed us to obtain multiple values of energy or energy loss for each particle combined with the 
CsI array. This enables us to identify the type of particle using the AE — E method. 

The Si-CsI detector is a versatile device that can be used for multi-alpha-particle detection. It combines the advantages of a 
Si detector, which can determine the position and energy of the particles with high precision, and a CsI scintillator. The SiPM 
is a key component of the detector, as it has high sensitivity and a fast response time, allowing for accurate detection and mea- 
surement of low-energy particles [43-45]. SiPMs coupled with scintillator crystals such as lutetium—yttrium oxyorthosilicate 
and LaBr3(Ce) crystals have been highly sought after in various fields [46, 47]. The combination of SiPMs and CsI exhibits 
favorable characteristics in terms of detection performance. The high sensitivity and small size of SiPMs enable detectors 
to achieve high spatial resolution and precise energy measurements. In addition, SiPMs have good single-photon resolution, 
providing high-precision photon counting and particle identification capabilities. 

One challenge in using DSSD is the possibility of multiple alpha particles simultaneously hitting the DSSD. The combination 
of the X position from the front side and Y position from the back side will lead to a multi-value for position measurement 
if more than one particle hits the DSSD, which makes it difficult to accurately calculate the momentum of the particle. To 
overcome this issue, an additional layer of DSSD can be used to obtain multi-channel signals and striking positions. Using a 
two-layer DSSD provides an additional set of positional information (other X and Y positions). The real hit position should 
be on the same line as the reaction center, which is the target. Based on this, we can determine the real positions of double-hit 
events. 

SiPMs are ideal for compact detectors owing to their high integration and small size. Previous studies have indicated that 
the spectral response of SiPMs is compatible with their CsI emission spectra. Because the SiPM+CsI can detect ~y, it can be 
used as a portable dosimeter for environmental dose equivalent and environmental dose rate equivalent measurements [45]. The 
combination used for proton detection consists of well-separated protons and deuterons at energies below 12 MeV. [48]. Using 
the Pulse Shape Discrimination (PSD) method, a particles can be easily separated from 6-y [51, 52]. In addition, the CsI crystal 
coupled to the SiPM was used to detect œ particles and low-energy protons and was tested using the Tandem Van-der-Graff of 
the INFN-Laboratori Nazionali del Sud (LNS) [49]. Thus, it can be inferred that the CsI-SiPM combination can also be applied 
in heavy-ion collision experiments. SiPMs are high-performance detectors that detect and amplify single-photon signals. The 
advantages of SiPMs over conventional photomultiplier tubes (PMTs) are manifold, such as high photon detection efficiency, 
single photon sensitivity, fast response time, low operating voltage, low power consumption, and small size. They also have 
good radiation and magnetic resistance, which make them adaptable to various environments and conditions. Therefore, using 
SiPMs as back-end photoelectric signal conversion devices for CsI can enhance detection capabilities by enabling efficient light 
signal readout and processing. 


I. SIMULATION 


Geant4 provides a toolkit to simulate the detection process, which offers various models for physical processes such as elastic 
scattering, inelastic scattering, ionization, and trajectory refinement. These models are essential for studying the interactions 
between particles and different materials in a detection design. Geant4 also allowed us to define the complex detector geometries 
in detail. This enabled us to accurately model the detector and evaluate its performance. 

The experimental results showed that the kinetic energy of a from the ?SSi breakup was lower than 300 MeV [10]. Alpha 
particles with energies below 520 MeV can be fully stopped in a 50 mm CsI crystal. Therefore, the energy of a particles can 
be measured using a CsI crystal. Furthermore, alpha-conjugate nuclei with higher masses, such as !7C, 160, ?°Ne, and 74Mg, 
are mostly stopped in DSSDs. Therefore, the experiment can detect the light and heavy particles of interest simultaneously. 


A. Angular Resolution 


The DSSDs can provide the position information of the particle hitting the detector, and the 0 and ¢ values of the particle 
emissions can be calculated using this position information. Assuming that the Z-axis is the direction of the incident particle 
beam, 0 and ¢ can be obtained as follows: 
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ọ = arctan( A (2) 


where x, y, z are the position information obtained from the DSSDs. 

Angle resolution is crucial for the performance of the system. To evaluate the desired angular resolution, we select several 
discrete values of 0 and ¢. The distance between the target and the detectors, which was 150 mm, was chosen to cover a forward 
angle of up to + 25°. The distribution of 0 without ¢ restriction is shown in Fig. 2(a). The detectors can cover a ¢ angle from 
0° to 360°, except for the gaps between the telescope arrays. The ¢-distribution is shown in Fig. 2(b). 
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Fig. 2. Distribution of a particles. Panels (a) and (b) show the @ and ¢ resolutions of the telescope array, which is placed 15 cm away from 
the target. Black curve represents the input value of the simulation, and the red and green curves represent the angular distribution calculated 
with single DSSD and double layer DSSDs, respectively. 


As Fig. 2 shows, the FHWM of angular resolution for 0 ranges from 0.2°to 0.4°, whereas the FHWM of angular resolution 
for ¢@ is approximately 2.2°. The measurements of 0 and ¢ show a noticeable difference in the angular resolution capability 
between the single DSSD and double-layer DSSD. The double-layer DSSD performs better in terms of angular resolution 
capability. These angular-resolution capabilities affect the reconstruction of fragmented nuclei during multiple alpha events. 
Hence, the use of two layers of DSSD enables more precise event reconstruction. The next section illustrates the differences in 
the reconstruction of the particle excited states between single- and double-layer DSSDs. 


B. Event Reconstruction 


By collecting and analyzing the data of the particles emitted from the nuclear collision process, we can infer the dynamics 
and nuclear state of the collision event. Event reconstruction is a crucial technique that allows physical information to be 
drawn from data. The performance and accuracy of event reconstruction are essential parameters affecting the performance of 
charged-particle telescope arrays. 

The outcome of low-energy particle collisions can be characterized by several parameters, such as the excitation energy of 
the fragment before it decays and the angular distribution of the fragment. These parameters provide valuable information on 
the dynamics and mechanisms of the collision process, as well as the properties and structures of the particles and fragments 
involved. Therefore, it is essential to measure and analyze the experimental parameters of low-energy particle collisions. 
Important and interesting physical parameters, such as Æ, denote the excitation energy of the projectile nuclei, and F,,; denotes 
the total kinetic energy of all fragments in the laboratory frame. 


and 
Etot = 5 Kiab + Krecoil (4) 


where Kem. denotes the kinetic energy in the center-of-mass frame of the summed particles. Q is the Q-value of a nuclear 
reaction, which represents the difference between the masses of the initial reactants and final products. Kj,» denotes the kinetic 
energy in the laboratory frame, and K,-¢co;; represents the kinetics of the recoil nuclear energy calculated by the conservation 
of energy and momentum. 

This information is crucial for gaining a more comprehensive understanding of the physical processes underlying an experi- 
ment. By combining the cascade decay process, we can gain insights into the fundamental properties of the particles involved 


as well as the energy distribution and dynamics of the system as a whole. In addition, according to the charge multiplicity dis- 
tribution, the temperature of the heavy-ion collisions can be determined, which can be applied to analyze the reaction dynamics 
[50]. 

This spectrum from ?°Si fragmented into 7 a particles [10] was used as an input for Geant4 simulation to model the detection 
process, and the alpha particles hit the detectors event-by-event. Such simulations are more realistic and allow for better 
evaluation and analysis of the performance of the detector under experimental conditions. However, based on these data, it is 
estimated that approximately 20% of the 7 alpha events may be lost because of the detection geometry efficiency. 
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Fig. 3. Panel (a) shows the excitation energy distribution of 7°Si reconstructed by 7 a particles. The panel (b) shows the distribution with a 
restriction on the input value between 83 MeV and 85 MeV. 


Fig. 3 shows the excitation energy distribution of 7°Si reconstruction from 7 a. It is evident that the calculated result with 
double DSSDs matches the simulated input. A single DSSD shows poor performance for the 7 a events. Because some events 
respond simultaneously on multiple silicon strips, this can interfere with the position resolution of a single DSSD, whereas 
double DSSDs can eliminate such effects through momentum analysis. 

This array of charged-particle detectors was specially designed and optimized for alpha particle detection. It can detect the 7 
q particle events that occurred in the experiment, as well as the fragmentation of alpha-conjugated nuclei into N alpha events. 
By analyzing the data from these events, we can reconstruct the excitation energies of the nuclei. 

One possible explanation for the occurrence of 7 a particles might be the cascade decay of 7°Si, which means that the decay 
product may contain other alpha-conjugated nuclei, such as Be and !?C. To investigate this hypothesis, the data from the 
alpha-conjugated nuclei should be analyzed. Reconstructing the cascade decay process of 7°Si to understand how it produces 7 
a particles is a novel approach for studying the nuclear structure and dynamics of 7°Si and its decay products. 
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Fig. 4. Excitation energy distribution of Be breaks up to 2 a particles and ‘°C breaks up to 3 a particles. 


The excitation energy of Be is low, and the 2 a particles from the decay have a very small emission angle in collisions, 
which requires a high angular resolution to detect the 2 a event. Be nuclei will be produced in large numbers in the reaction, 
indicating that there is a high rate of double-alpha hit events in the experiment. Therefore, it is important to simulate how the 


detector responds to the fragmentation events of Be. Fig. 4.(a) depicts the excitation energy spectra of 8 reconstructed using 
2 a-particles. When 2 a hit the detector simultaneously, two sets of X-Y values are obtained, corresponding to two sets of 
spatial positions. The single-layer DSSD cannot distinguish the real position of œ hit, which causes the result to deviate from 
the input value. Double-layer DSSDs provide another set of spatial positions from which the hit position can be determined 
through momentum analysis. In addition to Be, !?C can also be produced during the reaction. The production of !7C in 
stars via the triple-alpha process depends on a specific resonance of 1°C, called the Hoyle state. This resonance has important 
implications for nuclear structure theory and nuclear astrophysics applications. It is essential to evaluate the ability of the 
detectors to reconstruct the 1°C fragmentation events that occur in the Hoyle state. The simulation results are presented in 
Fig. 4(b). In the simulation for reconstructing !7C, the double-layer DSSDs outperformed the single-layer DSSD. As shown in 
Fig. 4, the charged-particle telescope array reconstructs the excitation energies of Be and !*C with reasonable accuracy. The 
alpha particles derived from the double-layer DSSD provided more accurate excitation energies than those derived from the 
single-layer DSSD. Therefore, a double-layer DSSD design is required for this purpose. 

Moreover, the results demonstrate that the telescope array is versatile and suitable for various experiments, in addition to 
those involving the disintegration of ?°Si into 7 a particles. It can also achieve good performance in experiments that generate 
multiple alpha particles or multiple fragments of different masses, owing to its high-angle resolution. 

The telescope array exhibited remarkable performance and accuracy in the excitation energy spectra of the reconstructed 
nuclei. Thus, it is a vital tool for studies involving these processes. 


IV. MEASUREMENT 


The simulation results indicate that the detector meets the performance requirements of the experiment; however, further 
validation is necessary to assess its actual performance. At present, the evaluation of the capability of the detector to detect 
charged particles is limited to the use of an 744Am source. The decay of 744Am primarily produces alpha particles with an 
energy of 5.486 MeV. These particles have a low penetration depth and can be blocked using a single-layer DSSD. Therefore, 
the energy resolution can be determined separately for the DSSD and CsI components. The energy resolution of the cesium- 
iodide scintillation detector was measured. 
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Fig. 5. Test results of 25 CsI+SiPM detectors with the 7*1Am a source. The data points are fitted with a Gaussian function. 


The energy spectrum of the telescope array consists of 25 CsI crystals coupled to 25 SiPMs with a 30 V voltage and gain of 
20 on the operation amplifier circuit, as shown in Fig. 5. It can be seen that there are small differences in performance for CsI 
crystals. This may be due to the variation in the alpha spectra because of the heterogeneity of the CsI crystals, although they 
are from one large block of CsI. In addition, the position of the radioactive source during the measurement may influence the 
results. 

Fig. 6 shows that the energy resolution is approximately 9%. Compared with the traditional design, sandwiched TiO2 can 
achieve comparable energy resolution [51, 52]. 

According to Fig. 7, the energy resolution of the CsI+SiPM composite detector depends on the voltage applied to SiPM. 
A higher voltage results in a higher gain, a smaller dynamic range, and better energy resolution. However, when detecting a 
particle with high energy, the increased noise and saturation effects caused by high voltage can affect the detector resolution. 
To detect high-energy, even up to 300 MeV alpha particles, the voltage applied to SiPM should be set as low as possible. 
However, the overall performance of the telescope array for high-energy particles, including the energy resolution, crosstalk, 
and charged-particle identification capability, must be assessed by performing online beam tests. 
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Fig. 6. Test results of 25 CsI+SiPM detectors. Channel number is the readout. 
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Fig. 7. Performance of the CsI+SiPM detector under different operating conditions. Channel numbers (red line) and energy resolution (black 
line) depend on the voltage applied to the SiPM. Channel number is the readout from MDPP-32. 


V. SUMMARY 


A novel detector array consisting of two DSSDs and a 5 x 5 cesium iodide scintillator array coupled with a SiPM array for 
particle identification and tracking of particle trajectories was designed and tested aiming at the detection of multiple a particle 
events with high position resolution. The high sensitivity and excellent timing performance of SiPM make it a good light signal 
readout module. CsI+SiPM detectors can be used in experiments to detect protons or deuterons with low energy. However, this 
is the first attempt to use the CsI+SiPM array to detect high-energy particles in heavy-ion collisions. The compact and robust 
telescope array is suitable for various experimental scenarios owing to the small size of SiPM and the sandwiching technology 
applied for CsI array. 


Using the Geant4 tookit, the detector responses to œ particles with different energies and angles have been simulated. The 
simulation results indicate that the detector array has high angular and energy resolutions, enabling accurate distinguishing 
of multiple a particle events. Additionally, the reconstruction result of ?°Si excitation energy using 7 a is satisfactory. This 
telescope array coupled with SiPM is a novel design for nuclear physics experiments, which can be used not only in experiments 
investigating exotic toroidal structures but also in future experiments to study the excitation states involving nuclear clusters, 
contributing to the exploration of the decay and dynamics of exotic nuclear shapes. 
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